Two Tautomers of Thiobarbituric Acid in One Crystal: The Experimental Charge Density Perspective

High-quality crystals of a certain polymorphic form of thiobarbituric acid containing both keto and enol tautomers in the asymmetric unit were obtained. High-resolution X-ray diffraction data up to sinθ/λ = 1.0 Å−1 were collected and subsequently successfully used for the refining of the multipolar model of electron density distribution. The use of a crystal containing both ketone and enol forms allowed a direct comparison of the topological analysis results and a closer look at the differences between these two forms. The similarities and differences between the deformation densities, electrostatic potentials, Laplacian maps and bond characteristics of the tautomers were analysed. Additionally, the spectrum of the intermolecular interactions was identified and studied from classical, relatively strong N-H···O and O-H···O hydrogen bonds through weaker N-H···S hydrogen bonds to weak interactions (for instance, C-H···O, C-H···S and N···O). The results of these studies point toward the importance of including both the geometrical features and the details of the electron density distribution in the analysis of such weak interactions.


Introduction
Tautomerism occurs when two or more constitutional isomers (with different connectivities) exist in a dynamical equilibrium. These forms can, therefore, (relatively) easily interconvert by transferring an atom (mainly hydrogen-prototrophic tautomerism) or even a group of atoms from one position of the molecule to another. This description does not introduce any clearly defined boundaries between isomers and tautomers. It is generally assumed that tautomers are isomers with a relatively low activation energy of interconversion of ca. 20 kcal/mol [1]. A certain type of this phenomenon, keto-enol tautomerism, seems to be the most appropriate for structural studies due to the slow transformation from one form to another and the possibility of the isolation of individual forms of tautomers [1]. The first example of tautomerism in a solid state was a case of keto-enol tautomers of acetyldibenzomethane isolated by Claaisen in 1896 [2]. A few structural studies of crystals with two stable tautomers in the unit cell have been published; for instance, the structures of isocytosine [3], L-His-Gly hemihydrate [4], N-(3-hydroxysalicylidene)-4-methoxy-aniline [5] and 4(5)-nitro-5(4)-methoxyimidazole [6].
It seems obvious that the use of a single sample containing two tautomers in one crystal would minimise the errors associated with the data collection, data processing and refinement and, therefore, would make the analysis of the generally fine differences between tautomers as efficient and detailed as possible. In the course of our ongoing studies of thioamines, we have realised that thiobarbituric acid (TBA) is a very good candidate for such an in-depth analysis of this phenomenon including a high-resolution electron density distribution analysis. It should be noted that the number of compounds containing two tautomers in an asymmetric unit of the unit cell of their crystal structure is extremely scarce and these structures often occur in the form of a solid solution, which effectively makes high-resolution research impossible.
There are two principal reasons favouring this choice: first, TBA can be obtained as well-defined, stable single crystals, meeting the criterion of high-resolution data collection and second, TBA is known to have a very rich (if not the richest) collection of "tautomeric polymorphs" [7] including one crystal structure of a pure enol form, three crystal structures of keto isomers and one co-crystal of a keto/enol form, which tend to be the thermodynamically most stable form (at room temperature) and include the highest number of hydrogen bonds in the crystal structure [7]. Additionally, thiobarbituric acid is the biological relevant compound, often used in measuring the lipid peroxidation products in cells, tissues, body fluids and food [8,9].

X-ray Diffraction
A TBA polymorph containing both the keto and enol forms of thiobarbituric acid was crystallised from an acetonitrile solution by very slow evaporation. Regular pink crystals suitable for X-ray data collection were obtained. The diffraction data were collected by the ω-scan technique using graphite-monochromated MoK α radiation (λ = 0.71073 Å) at 100 (1) K on a Rigaku XCalibur four-circle diffractometer with an EOS CCD detector. The data were corrected for Lorentz polarisation as well as for absorption effects [10]. Precise unit cell parameters were determined by a least-squares fit of the reflections of the highest intensity chosen from the whole experiment. The structure was solved with SHELXT [11] and refined (at an independent atom model approximation) with the full matrix leastsquares procedure on F 2 by SHELXL [12] within the Olex-2 suite of programs [13]. All non-hydrogen atoms were anisotropically refined. Hydrogen atoms were found in the difference Fourier maps and freely refined. The crystal data (Table 1) showed that it was the polymorph described by Chierotti et al. [7] as polymorph IV: it crystallised in a P2 1 /c space group with two molecules in an asymmetric unit (keto and enol form) and the unit cell parameters were similar; the differences could be related to the different experiment temperatures (room temperature in [7], 100 K in these studies). The thermalellipsoid representations of both molecules are presented in Figure 1 and the most important crystallographic parameters are summarised in Table 1. ing two tautomers in an asymmetric unit of the unit cell of their crystal structure is extremely scarce and these structures often occur in the form of a solid solution, which effectively makes high-resolution research impossible. There are two principal reasons favouring this choice: first, TBA can be obtained as well-defined, stable single crystals, meeting the criterion of high-resolution data collection and second, TBA is known to have a very rich (if not the richest) collection of "tautomeric polymorphs" [7] including one crystal structure of a pure enol form, three crystal structures of keto isomers and one co-crystal of a keto/enol form, which tend to be the thermodynamically most stable form (at room temperature) and include the highest number of hydrogen bonds in the crystal structure [7]. Additionally, thiobarbituric acid is the biological relevant compound, often used in measuring the lipid peroxidation products in cells, tissues, body fluids and food [8,9].

X-ray Diffraction
A TBA polymorph containing both the keto and enol forms of thiobarbituric acid was crystallised from an acetonitrile solution by very slow evaporation. Regular pink crystals suitable for X-ray data collection were obtained. The diffraction data were collected by the ω-scan technique using graphite-monochromated MoKα radiation (λ = 0.71073 Å) at 100 (1) K on a Rigaku XCalibur four-circle diffractometer with an EOS CCD detector. The data were corrected for Lorentz polarisation as well as for absorption effects [10]. Precise unit cell parameters were determined by a least-squares fit of the reflections of the highest intensity chosen from the whole experiment. The structure was solved with SHELXT [11] and refined (at an independent atom model approximation) with the full matrix leastsquares procedure on F 2 by SHELXL [12] within the Olex-2 suite of programs [13]. All non-hydrogen atoms were anisotropically refined. Hydrogen atoms were found in the difference Fourier maps and freely refined. The crystal data (Table 1) showed that it was the polymorph described by Chierotti et al. [7] as polymorph IV: it crystallised in a P21/c space group with two molecules in an asymmetric unit (keto and enol form) and the unit cell parameters were similar; the differences could be related to the different experiment temperatures (room temperature in [7], 100 K in these studies). The thermal-ellipsoid representations of both molecules are presented in Figure 1 and the most important crystallographic parameters are summarised in Table 1.

Multipolar Modelling
The multipolar model of the electron density distribution was refined against structure factor amplitudes with MoProSuite [14] using the Hansen-Coppens model [15], in which the pseudoatom electron density is described as: ρ atom (r) = ρ core (r) + P val κ 3 ρ val (κr) + Σ l κ' 3 R l (κ'r)Σ m P lm± (θ, ϕ) where P val denotes the valence population, P lm± are the multipole populations and κ and κ' are the contraction/expansion parameters, respectively. R l denotes the radial Slater-type function. The oxygen, carbon and nitrogen atoms were refined up to the octupolar level (l max = 3), the sulphur atoms up to the hexadecapole level (l max = 4) and the hydrogen atoms up to the dipole level (l max = 1). The n l and ζ l values were respectively equal to 4, 4, 4, 6, 8 and 3.851 au −1 for the S atoms; 2, 2, 2, 3 and 4.47 au −1 for the oxygen atoms; 2, 2, 2, 3 and 3.84 au −1 for the nitrogen atoms; and 2, 2, 2, 3 and 3.18 au −1 for the carbon atoms. In the case of the hydrogen atoms, the n l and ζ l values were 1, 1 and 2.00 au −1 . The core and valence scattering factors were calculated from Clementi wave functions [16]. The reflections up to sinθ/λ = 1.0 Å −1 were used in the refinement. The ADPs for the H atoms were estimated using the SHADE server [17]. The X-H distances were initially constrained to the average values obtained from neutron diffraction studies and refined at the final stages of the multipolar refinement. The third-order anharmonic model of thermal motion was implemented for the S2 and S12 atoms. An I/σ cut-off = 2.0 was applied for the final refinement.

Verification of the Model Correctness
A residual density analysis (RDA) [18] was performed for the final models. The plots representing the results of this analysis are included in the Supplementary Materials ( Figure S1) and they show that, in the case of both resolutions, the parabolas were symmetrical and regular, indicating the appropriate quality of the final models, which was also confirmed by the clear residual electron density maps (included in the Supplementary Materials). Additionally, the Hirshfeld rigid bond test [19] was applied to the covalent bonds after the final refinement. For each pair of non-hydrogen atoms, the values of ∆Z AB 2 (differences between the components of the Uij tensor along the bond) did not exceed the acceptable value of 10 −3 Å 2 . The largest differences were found for N13-C14 (0.00048) and for C15-C16 (0.00045). The tables containing all values for the rigid bond test as well as the residual electron density maps and dynamic deformation electron density maps are included in the Supplementary Materials (Table S1, Figures S2 and S3).

Electron Density Maps
The comparison of the (multipolar) electron density distribution in the two tautomers started from the final static deformation 2D and 3D electron density maps ( Figure 2). One would expect the differences in the electron density distribution in the heterocyclic rings as the result of a significant electronic change. This really was the case: in the keto tautomer, the electron density distribution within the ring was less symmetrical due to a loss of the full aromaticity caused by the sp 3 hybridisation of the C15 atom. This atom was no longer in the plane of the ring and the deviation of this atom from the least-squares plane calculated for the remaining five ring atom was equal to 1.171 (3) Å. In the keto tautomer, the distribution of the electron density around the sulphur (S12) and oxygen (O14) atoms was also less symmetrical than in the enol form. after the final refinement. For each pair of non-hydrogen atoms, the values of ΔZAB 2 (differences between the components of the Uij tensor along the bond) did not exceed the acceptable value of 10 −3 Å 2 . The largest differences were found for N13-C14 (0.00048) and for C15-C16 (0.00045). The tables containing all values for the rigid bond test as well as the residual electron density maps and dynamic deformation electron density maps are included in the Supplementary Materials (Table S1, Figures S2 and S3).

Electron Density Maps
The comparison of the (multipolar) electron density distribution in the two tautomers started from the final static deformation 2D and 3D electron density maps ( Figure 2). One would expect the differences in the electron density distribution in the heterocyclic rings as the result of a significant electronic change. This really was the case: in the keto tautomer, the electron density distribution within the ring was less symmetrical due to a loss of the full aromaticity caused by the sp 3 hybridisation of the C15 atom. This atom was no longer in the plane of the ring and the deviation of this atom from the least-squares plane calculated for the remaining five ring atom was equal to 1.171 (3) Å. In the keto tautomer, the distribution of the electron density around the sulphur (S12) and oxygen (O14) atoms was also less symmetrical than in the enol form.

Bond Critical Points (Intramolecular)
With a high-quality electron density distribution, it was possible to apply the atoms in molecules approach [20] in order to calculate the critical points of this distribution; in

Bond Critical Points (Intramolecular)
With a high-quality electron density distribution, it was possible to apply the atoms in molecules approach [20] in order to calculate the critical points of this distribution; in particular, the bond critical points (BCP) of the signature (3, −1). These data (Table 2, Figure 3) allowed for a more detailed analysis of the differences between the tautomers at the level of intramolecular electron density distribution. The full characteristics of the critical points are presented in the Supplementary Material. particular, the bond critical points (BCP) of the signature (3, −1). These data ( Table 2  The highest value of the electron density in the BCP was observed (as expected) for the C=O bond in both tautomers but the appropriate value for the single C6-O6 bond in the enol tautomer was significantly smaller (Table 2). This was accompanied by differences in the neighbouring C-C bond; the value of the electron density at the CP of the C5-C6 bond (enol) was much higher than that of the C15-C16 bond (keto). The value for C5-C6 was even larger than several C-N bonds; this tendency expanded toward the next bond, C4-C5. Another expected change was related to the loss of the ring aromacity in the keto tautomer and the influence of this loss on the delocalisation of the electrons of the C=O and C=S bonds. As expected, larger values of electron density at these bonds were observed in the keto tautomer as an effect of the above-mentioned change. As a further consequence, a noticeable decrease of the electron density on the C12-N13 and C12-N11 bonds (keto) in comparison with the C2-N3 and C2-N1 (enol) bonds was also observed.  The highest value of the electron density in the BCP was observed (as expected) for the C=O bond in both tautomers but the appropriate value for the single C6-O6 bond in the enol tautomer was significantly smaller (Table 2). This was accompanied by differences in the neighbouring C-C bond; the value of the electron density at the CP of the C5-C6 bond (enol) was much higher than that of the C15-C16 bond (keto). The value for C5-C6 was even larger than several C-N bonds; this tendency expanded toward the next bond, C4-C5. Another expected change was related to the loss of the ring aromacity in the keto tautomer and the influence of this loss on the delocalisation of the electrons of the C=O and C=S bonds. As expected, larger values of electron density at these bonds were observed in the keto tautomer as an effect of the above-mentioned change. As a further consequence, a noticeable decrease of the electron density on the C12-N13 and C12-N11 bonds (keto) in comparison with the C2-N3 and C2-N1 (enol) bonds was also observed.

Atomic Charges
The atomic charge distribution, calculated here by the integration of the electron density distribution over the atomic basins and defined by zero-flux surfaces, was generally quite similar in both tautomers ( Table 3). The main differences, corresponding with the differences in the electron densities at the bond critical points described above, were observed for the pair C6/C16; the carbon atom from the keto tautomer was more positive than the corresponding one from the enol form, which could be rationalised by the observation that the oxygen atom in the keto form (C=O) much better stabilised the negative charge than the atom in the hydroxyl group. On the other hand, it should be noted that the differences between the O6 (enol) and O16 (keto) atoms were more significant than between the carbon atoms. However, for this pair, one should compare the whole hydroxyl group O6H6 with O16 and not just the oxygen atoms. Once this was done, the hydroxyl group still had a negative charge but it was definitely smaller (0.325 e) than for O16 atom 1.061 e. A similar reasoning could be made for the C15 (keto, sp 3 ) and C5 (enol, sp 2 ) atoms, which taken alone had the second highest values of charge difference but after adding the bonded hydrogen atoms, both groups had similar charges of 0.11 and 0.12 e.

Electrostatic Potential
The electrostatic potential distribution for both tautomeric forms was basically similar (Figure 4) although the enol tautomer was definitely more electrophilic in the ring area. In both cases, the ketone oxygen atoms defined the most electronegative parts of the molecules and the sulphur atoms were rather neutral. The differences concerned mainly the ring regions with an almost exclusively positive potential in the case of the enol tautomer and, as a more complicated picture, divided into zones of a positive and neutral potential for the keto tautomer.

Intermolecular Interactions
It was found that thiobarbituric acid had a very rich collection of polymorphic forms; this was related to the presence and availability of very good donors and acceptors of hydrogen bonds in the molecule allowing for the constructions of different stable crystal architectures. It was also found that the most stable form at room temperature was the

Intermolecular Interactions
It was found that thiobarbituric acid had a very rich collection of polymorphic forms; this was related to the presence and availability of very good donors and acceptors of hydrogen bonds in the molecule allowing for the constructions of different stable crystal architectures. It was also found that the most stable form at room temperature was the form with two tautomers in the unit cell [7], the one studied here. With very good data and a more advanced model of the electron density distribution, we could attempt to analyse in detail the factors that at least could be responsible for this stability. As many as 22 intermolecular bond critical points were found in this structure; the selected characteristics of these points are listed in Table 4 and the full collection of data for the intermolecular BCP can be found in the Supplementary Materials (Table S3).
There were nine interactions that met the criterion of the presence of a hydrogen bond, as defined by Munshi and Guru Row [21] and three of them could be classified as strong HBs (N-H···O and O-H···O). The distinct second-in-strength group contained two N-H···S hydrogen bonds and then there was a continuum of weaker and weaker interactions (C-H···O, C-H···S, etc.). Only one of the contacts, C-H···S (namely, cp 45), seemed to correspond more with the van der Waals-type interaction (region 3) than to the hydrogen bonding but even in this case, the interaction seemed to be important and led to the occurrence of the appropriate ring critical point.
Two of the strongest hydrogen bonds (both in terms of their topological properties and of their geometries (Table 3)), together with weaker but well-defined N-H···S and C-H···S interactions, formed layers of alternating tautomers ( Figure 5). tions (C-H•••O, C-H•••S, etc.). Only one of the contacts, C-H•••S (namely, cp 45), seemed to correspond more with the van der Waals-type interaction (region 3) than to the hydrogen bonding but even in this case, the interaction seemed to be important and led to the occurrence of the appropriate ring critical point.
Two of the strongest hydrogen bonds (both in terms of their topological properties and of their geometries (Table 3)), together with weaker but well-defined N-H•••S and C-H•••S interactions, formed layers of alternating tautomers ( Figure 5). It should be noted, however, that the geometrical parameters (Table 5) did not explain the diversity of the electron density and Laplacian values at the BCPs corresponding with the N-H ••S hydrogen bonds as both the donor acceptor distances and the D-H•••A angles were very similar. The analysis of the appropriate maps of the Laplacian and electron density distribution indicated that these differences could be connected to the much higher charge concentration in the direction of the hydrogen bond in the case of the S12 atom (keto form) than for the S2 (enol), which in turn resulted in a stronger N-H•••S interaction ( Figure 6). It should be noted, however, that the geometrical parameters (Table 5) did not explain the diversity of the electron density and Laplacian values at the BCPs corresponding with the N-H ··S hydrogen bonds as both the donor acceptor distances and the D-H···A angles were very similar. The analysis of the appropriate maps of the Laplacian and electron density distribution indicated that these differences could be connected to the much higher charge concentration in the direction of the hydrogen bond in the case of the S12 atom (keto form) than for the S2 (enol), which in turn resulted in a stronger N-H···S interaction ( Figure 6).    In addition to these interactions, which were definitely the most important for the determination of the crystal architecture, there were also a number of weak or even very weak but still interesting interactions between the thioamide sulphur atoms and the delo-calised electrons of the thioamide groups (cp34-37). The importance of these interactions (or the geometrical necessity of their existence) seemed to be confirmed by the fact that the sulphur atoms from both tautomeric forms were involved in such contacts (Tables 3 and 4; Figure 7). C5-H5 S12 1.08 2.82 148 3.7822 (6) x, −y + 1/2, z + 1/2 N13-H13 S2 1.00 2.31 164 3.27662(9) − x + 1, −y, −z C15-H15A S2 1.09 2.97 136 3.8254 (6) x − 1, −y + 1/2, z + 1/2 In addition to these interactions, which were definitely the most important for the determination of the crystal architecture, there were also a number of weak or even very weak but still interesting interactions between the thioamide sulphur atoms and the delocalised electrons of the thioamide groups (cp34-37). The importance of these interactions (or the geometrical necessity of their existence) seemed to be confirmed by the fact that the sulphur atoms from both tautomeric forms were involved in such contacts (Tables 3  and 4; Figure 7). In detail, the sulphur atom S12 (keto tautomer) interacted with the carbon atom C12 (1 − x, −y, −z) and the same sulphur atom was involved in the S•••N (1 − x, −½ + y, −½ − z) contact (cp37). In turn, the sulphur atom S2 (enol) was also involved in similar contacts but in this case the electron density around the S2 atom, which could be ascribed to the In detail, the sulphur atom S12 (keto tautomer) interacted with the carbon atom C12 (1 − x, −y, −z) and the same sulphur atom was involved in the S···N (1 − x, − 1 2 + y, − 1 2 − z) contact (cp37). In turn, the sulphur atom S2 (enol) was also involved in similar contacts but in this case the electron density around the S2 atom, which could be ascribed to the free electron pairs, was much more deflected in the direction of these interactions. This detail, however, did not significantly affect the parameters characterising the appropriate critical points, which were comparable for both thiobarbituric acid molecules involved in the interactions.
The other interesting contact, for which the interaction path and a critical point could be determined, was found between the nitrogen atoms of the heteroatomic ring and the oxygen atoms of both tautomeric forms (Figure 8). These interactions were homomolecular (i.e., between the same tautomeric forms) but were found for both forms; there were ketoketo and enol-enol pairs in the structure. In both cases, the mutual arrangement of the rings was antiparallel with a significant shift of molecules.
Materials 2022, 15, x FOR PEER REVIEW 10 of 12 free electron pairs, was much more deflected in the direction of these interactions. This detail, however, did not significantly affect the parameters characterising the appropriate critical points, which were comparable for both thiobarbituric acid molecules involved in the interactions. The other interesting contact, for which the interaction path and a critical point could be determined, was found between the nitrogen atoms of the heteroatomic ring and the oxygen atoms of both tautomeric forms (Figure 8). These interactions were homomolecular (i.e., between the same tautomeric forms) but were found for both forms; there were keto-keto and enol-enol pairs in the structure. In both cases, the mutual arrangement of the rings was antiparallel with a significant shift of molecules.  Interestingly-and to an extent, strangely-in this interaction, another (3, −1) bond critical point together with the corresponding bond paths were found in the locations where one would expect ring critical points (3, + 1) (cp48, cp46 (Figure 8)). These BCPs were related to the N•••O interactions and had quite high values of electron density as well as Laplacian and energy density at the critical point so their influence on the crystal packing could be quite significant. It should be stressed that the pure analysis of only the geometrical features (distances, angles) would not give hints for considering this type of contact as relevant; however, the detailed analysis of the 3D electron density maps led to the assumption that this interaction had an electrostatic nature due to the fact that the areas of charge concentration were paired up with the neighbouring counterparts of electron depletion. It could be stated that these observations actually suggested that the distance  Interestingly-and to an extent, strangely-in this interaction, another (3, −1) bond critical point together with the corresponding bond paths were found in the locations where one would expect ring critical points (3, + 1) (cp48, cp46 (Figure 8)). These BCPs were related to the N···O interactions and had quite high values of electron density as well as Laplacian and energy density at the critical point so their influence on the crystal packing could be quite significant. It should be stressed that the pure analysis of only the geometrical features (distances, angles) would not give hints for considering this type of contact as relevant; however, the detailed analysis of the 3D electron density maps led to the assumption that this interaction had an electrostatic nature due to the fact that the areas of charge concentration were paired up with the neighbouring counterparts of electron depletion. It could be stated that these observations actually suggested that the distance between the rings was not so crucial for the importance of the contact as their mutual arrangement (Figure 9). Despite the similar distance of the rings, the critical points were found only in the case of the antiparallel orientation of the rings.

Conclusions
An in-depth analysis of the experimentally determined electron density in crystals containing two different tautomers of thiobarbituric acid was performed. The electron density distribution was obtained by means of the multipolar modelling of the molecular model using high-resolution X-ray diffraction data. The differences between the keto and enol tautomers were clearly seen in the deformation density maps as well as in the Laplacian and electrostatic potential maps. The topological analysis (atoms in molecules) of the gradient of the electron density allowed for the determination of the fine differences between the bond critical points in both tautomers. Additionally, this analysis allowed also for a detailed study of the intermolecular interactions, starting from well-defined strong hydrogen bonds through more obscure N-H•••S bonds to weak interactions; for instance, the C-H•••S, C-H•••O or N•••O types. It should be stressed that the results of these investigations showed that for the intermolecular interactions one should consider not only the geometrical parameters of the contacts but also the details-sometimes quite fine-of the electron density distribution. Only such a combined insight (maybe enriched by the energy calculations) could give the whole picture of the studied interaction. Separate analyses of the geometry or electron density distribution resulted in an under-or overestimation of the importance of a particular contact.

Conclusions
An in-depth analysis of the experimentally determined electron density in crystals containing two different tautomers of thiobarbituric acid was performed. The electron density distribution was obtained by means of the multipolar modelling of the molecular model using high-resolution X-ray diffraction data. The differences between the keto and enol tautomers were clearly seen in the deformation density maps as well as in the Laplacian and electrostatic potential maps. The topological analysis (atoms in molecules) of the gradient of the electron density allowed for the determination of the fine differences between the bond critical points in both tautomers. Additionally, this analysis allowed also for a detailed study of the intermolecular interactions, starting from well-defined strong hydrogen bonds through more obscure N-H···S bonds to weak interactions; for instance, the C-H···S, C-H···O or N···O types. It should be stressed that the results of these investigations showed that for the intermolecular interactions one should consider not only the geometrical parameters of the contacts but also the details-sometimes quite fine-of the electron density distribution. Only such a combined insight (maybe enriched by the energy calculations) could give the whole picture of the studied interaction. Separate analyses of the geometry or electron density distribution resulted in an under-or overestimation of the importance of a particular contact.
Supplementary Materials: The following supporting information can be downloaded at https: //www.mdpi.com/article/10.3390/ma15010223/s1, Figure S1: RDA plot, Figure S2: Dynamic deformation density maps, Figure S3: Residual electron density maps, Table S1: Rigid bond test, Table S2: Bond critical points, Table S3: Bond critical points for intermolecular interactions. Funding: These studies were supported by the "Scholarship support for doctoral students in fields considered strategic from the point of view of Wielkopolska's development".